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Abstract 

We study the implication of finite squark and gluino decay widths for threshold 
resummation of squark and gluino production cross sections at the LHC. We find 
that for a moderate decay width (T/fh < 5%) higher-order soft and Coulomb 
corrections are appropriately described by NLL calculations in the zero-width 
limit including the contribution from bound-state resonances below threshold, 
with the remaining uncertainties due to finite-width effects of a similar order as 
the ambiguities of threshold-resummed higher-order calculations. 



1 Introduction 



A recent development in the effort to improve theoretical predictions for the pair production 
of squarks and gluinos at hadron colliders has been the calculation of higher-order threshold 
corrections in QCD [l}jl~2], i.e. corrections that are enhanced in the limit of a small relative 
velocity (3 of the produced squark/gluino pair 







Am? 



->0, 



with the partonic centre-of-mass energy s and the average mass m of the sparticle pair. 
These threshold corrections are given by two contributions, Coulomb corrections, that 



occur in powers of a s j /3, and threshold logarithms of the form a™ ln m f3. In the limit ( 1.1 ) 



the expansion in the strong coupling constant must be reorganized, since both types of 
corrections can become of order one, 



a s In j3 
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Resumming all enhanced contributions to all orders using eikonal methods 13 16 and 
non-relativistic field theory (TtJ , the cross section assumes the form 



r (0) 
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exp 



In f3 g (a s In f3) + gi(a s In (3) + a s g 2 {a s In /?) + ... 
(NLL) 
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For squark and gluino production, the corrections beyond next-to-leading order (NLO) can 
be large, becoming as large as the NLO cross sections for the case of heavy gluinos. 

In the calculations of threshold corrections the produced squarks and gluinos are usu- 
ally treated as stable, while in most realistic SUSY models they are highly unstable and 
decay promptly]]] In a more realistic treatment including finite decay widths T, the thresh- 
old singularities are smoothed [17] since the decay width sets the smallest perturbative 
scale in the process. This raises the question of whether the resummation based on the 



counting (1.2) is appropriate in the presence of the new scale T and to which extent the 



large higher-order threshold corrections found for stable squarks and gluinos provide a 
reasonable approximation to the more realistic unstable case. One effect of the finite de- 
cay width is the emergence of a non-vanishing partonic cross section below the nominal 
production threshold s = 4m 2 , as discussed for stop-pair production in 20 , and more re- 



cently for gluino pair 21 22 and squark-gluino production 23 . In our previous combined 



resummation of threshold logarithms and Coulomb corrections for squark and gluino pair- 



production 11,12 the contribution from the region below the nominal threshold to the 



1 We do not consider the case of stops and gluinos that are stable on collider timescales. For recent 



work on higher-order corrections in this case see 18 19 . 
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total production cross section has been taken into account in the zero- width limit, where 
the smooth invariant-mass distribution turns into a series of would-be bound-state poles. 

As we will show in this note, for moderate decay widths T/rh < 5%, the finite-width 
effects on higher-order QCD corrections are small and within the error estimate of the NLL 
cross sections with soft and Coulomb resummation and bound-state corrections. This cov- 
ers the relevant production channels in the MSSM, since larger decay widths only arise for 
heavier gluinos with a large gluino-squark mass splitting, in which case gluino production 
is suppressed compared to squark production. Therefore the treatment of NLL soft and 
Coulomb corrections in [12] captures the dominant higher-order QCD effects also when the 
instability of squarks and gluinos is taken into account. Furthermore, the dominant effect 
of the finite decay widths arises at NLO, so the impact on higher-order corrections is small. 

It should be emphasized that the aim of our study is the assessment of finite-width 
effects on QCD corrections beyond LO rather than a complete treatment at leading and 
next-to-leading order. Studies of other aspects of finite squark and gluino decay widths 



include that of non-resonant tree diagrams 24 ; the accuracy of the narrow-width approxi- 



mation for small mass-splittings of the decaying particle and its decay products 25 26 , the 



consistent treatment of off-shell effects in Monte-Carlo programs 27 and NLO corrections 



to production and decay for the squark-squark production process 28 



This article is organized as follows. In Section [2] we specify the considered production 



and decay processes and review the effective-theory treatment of unstable particles 29 - 32 



and soft and Coulomb resummation near threshold [11]. In Section |3| we estimate the 
impact of the screening of threshold corrections by the finite decay width as well as the 
magnitude of non- resonant corrections. In Section [4] we study the invariant mass spectrum 
for squark and gluino production processes at the LHC and the total production cross 
sections. We find that for moderate decay widths the finite-width corrections to the total 



cross section are within the residual uncertainty of the NLL calculation of 12 



2 Theoretical framework 

2.1 Production and decay of squarks and gluinos in SQCD 

In this note we consider the pair-production processes of squarks q and gluinos g at hadron 
colliders in supersymmetric QCD (SQCD) that proceed through the partonic production 
channels 

pp'^hhX, (2.1) 

where Sj G {q, q, g} denote the two produced sparticles and p,p' G {q, q, g} the initial-state 
partons. The average mass of the sparticle pair is given by 

m= m5l + m *. (2.2) 

In general, the produced squarks and gluinos can decay via long cascade decay chains to a 
final state of standard model particles and the lightest supersymmetric particle (assuming 
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the latter is stable on detector time scales). For this initial study of finite- width effects we 
limit ourselves to squark and gluino decay in pure SQCD, i.e. we neglect the decay width 
of the lightest coloured sparticle and consider the decay modes 

g -> qg, q -> qg, rriq > m s , 

9 -> qq, 9 -> qq, m ? - < m § . (2.3) 

We treat the masses of the different squark flavours as degenerate and equal to rriq, while 
the gluino mass equals m g ^ rriq. The full one-loop SUSY-QCD partial widths of the 
above decay processes have been calculated in (33]. The NLO corrections to the LO partial 
widths were found to be of the order of 30 — 50% and will not be given here. In order to 
estimate the finite-width effects we will employ the LO tree level SUSY-QCD widths, 



r ( ^ 9 §) = ^i(i-(!=4) 

2 V \m„-/ 



2\ 2 

, rriq > m- g , 



vq, 

r(9 -+ qq,qq) = °^f^(l - (^fT^ < ™f ( 2 -4) 
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Since the numerical values are only used for illustration of the finite-width effects, we 
only consider gluino decay into light-flavour and bottom squarks that are treated as mass- 
degenerate, i.e. we set nf — 5. Therefore we neglect the decay into stops and the related 
additional dependence on the parameters of the stop sector unless stated otherwise. The 



effect of gluino decay to stops on our results is discussed briefly in Section |4.2| In the full 
MSSM, the lighter coloured sparticles are, of course, unstable with respect to electroweak 
two- or three-body decays q — > xq an d 9 qqx f° r charginos or neutralinos Xi which have 



been calculated in leading order in 34 , 35 and for which state-of-the-art predictions are 



implemented in 36 -39 . These effects would be straightforward to include in our framework 



but would introduce a dependence on the electroweak parameters of the MSSM. Therefore 



we will use the tree-level SQCD decay widths (2.4) for the purpose of estimating the size of 
finite-width effects in higher-order QCD corrections. Since the electroweak decay widths 
are typically of the order Y/fh < 1%, this should not affect our qualitative conclusions. 
We therefore consider production and decay processes of squarks and gluinos at hadron 
colliders of the type 

pp' -> 5iS 2 -»■ (s' 2 p 2 ), (2.5) 

with massless partons pi/2- Such four-body final states arise in SQCD for gluino pair 
production for the case m g > rriq or squark- ( ant i-) squark production for rriq > rrig. The 
treatment of a three-body final state that arises for squark-gluino production in our ap- 
proximation should be obvious. In detail, the partonic production and decay processes of 
squarks and gluinos for the two mass hierarchies are as follows: 

m-g < rriq: 99,qiqi^99, 

qi9 ^qg -> qgg, 

gg, qiCjj -> qq -> qqgg, 
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nig > rriq 



QiQj ^qq^ qqgg, 

99, q^j -> qq, 
q^j qq, 

qi9 ^qg -> gg?, 
00,9** -> qq ^ qqqq 



(2.6) 



where z, j = u, d, s, c, b and all charge conjugated processes are understood. 

The complete gauge-invariant set of Feynman diagrams for the production of the three- 



or four-body final states of "stable" particles (2.6) contain doubly-resonant diagrams, i.e. 



diagrams where the particles indicated in the first step of (2.6) appear as internal lines, and 



singly- or non-resonant diagrams, where only one or none of the unstable squarks or gluinos 
is present. Examples of these topologies are given in Figure [TJ Beyond leading order, the 
production processes cannot be strictly separated since, e.g. the process qig — > qqq enters 
the real NLO corrections to gg-pro duct ion. Furthermore, the non-resonant diagrams can 
contain collinear singularities (e.g. the last topology given in Figure [T]) that would cancel 
against higher-order virtual corrections and require taking kinematic cuts on the final state 
particles into account. 



A full NLO QCD calculation of a process comparable to (2.5) including finite-width 
effects has been performed in the standard model for the case of bbW + W~ production 40 



41 , but not yet for an MSSM process. Very recently, a tool for automatic calculation of 

, but so far applied to the on-shell 



MSSM processes at NLO has been presented in 



processes (2.1). Short of a full NLO calculation, a well-defined approximation method for 



calculating radiative corrections to production and decay processes is based on the pole 
In this framework the total partonic production cross sections for the 



expansion 



processes (2.6 



43 



44 



is written in the form 



a 



i "non-res / ~\ 
+ <7pp> ( S J 



(2.7) 



The doubly-resonant cross-section a res is defined in a consistent way by an expansion of 
the S'-matrix around the complex poles of the Dyson-resummed propagators of the un- 
stable particles 43 , 44 . At NLO it contains both factorizable corrections to production 



and decay, as well as non-factorizable corrections connecting production, propagation and 
decay stages. Non-factorizable corrections related to the final state cancel for the total 
cross section 45,46 , but not for differential distributions. The fully differential factoriz- 



able corrections to production and decay have recently been computed for squark-squark 
production [28]. For the remaining processes, at present only the corrections to the to- 
tal production cross sections 47 and decay rates 33 are available. The "non-resonant" 
cross-section <j non - res includes singly- and non-resonant Feynman-diagrams as well as con- 
tributions from doubly-resonant diagrams where one or both of the si or §2 lines are far 
off-shell. 
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Figure 1: Examples of Feynman-diagram topologies contributing to the generic production 
and decay process (2.5). The second non-resonant diagram contains collinear singularities. 



2.2 Effective theory framework for unstable particles 



For the investigation of the interplay of higher-order threshold corrections in QCD and 
finite squark and gluino lifetimes, we are interested in the partonic cross sections for the 
processes (2.6) near the partonic production threshold (1.1). For the precise definition of 



the resonant and non-resonant contributions to the cross section we adopt the effective- 
theory approach to unstable particles 29 30 , generalizing the treatment of 1^-boson pair 



production at an electron positron collider |31[|32| . The aim of the effective-theory approach 
is to provide a precise prediction of the partonic cross section for partonic centre-of-mass 
energies in the vicinity of the threshold, s-4m 2 ~ mT. This is achieved by a simultaneous 
expansion in the quantity 

s — 4rh 2 



m? 



and the coupling constants. For power-counting purposes, we therefore set 

13 ~ {T/rh) 1 ' 2 ~ 8 1 ' 2 . 



2.8 



(2.9) 



The leading-order resonant production cross-section for .S-wave production processes is of 
the form 



<7pp/(S 



a 



(2.10) 



Near the production threshold, the processes (2.5) can be treated in an effective theory 
where the light-partons are described by collinear fields <ft, <ft' in soft-collinear effective 
theory (SCET) (48- 



50 



and the pair-produced squarks and gluinos are described by non- 
relativistic fields ip,ij)' in potential non-relativistic QCD (pNRQCD) [51] . In the leading 
effective Lagrangian, the collinear and non-relativistic fields interact only through the 
exchange of soft gluons, which give rise to the non-factorizable corrections mentioned above. 
This effective theory has also been used for soft-gluon and Coulomb resummation in 
where more details can be found. 
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In the EFT framework, the total partonic cross section for the process (2.5 ) is computed 
from the imaginary part of the partonic forward-scattering amplitudes A pp > of the processes 



5 



iA 



pp 




Figure 2: Diagrammatic representation of the leading resonant and non-resonant contri- 
bution to the forward-scattering amplitude in unstable particle effective theory. 




+ 




Figure 3: Computation of the matching coefficient of the production operator for squark- 
antisquark production from the quark-antiquark initial state. 





Figure 4: Computation of the matching coefficient of the four-parton operator for squark- 
antisquark production and decay. The first diagram is an example of a double-resonant 
diagram, the second of the interference of a double-resonant and a single-resonant diagram. 



pp' — > pp' which read^] 



d A x (pp'\T [iOl od (0)iO piod (x)] \pp') + (pp'\tO 4p (0)\pp'} . (2.11) 

The "production operator" P rod = Cprod 4>4>' V* V describes the resonant production of 
the 5iS2-pair from pp' collisions while the expectation value of the four-parton operators 
describes the non-resonant contributions. The two terms in (2.11) 



O 



c 



ip 



therefore correspond to the resonant and non-resonant cross sections defined earlier in (2.7) 



see Figure [2] The cross section for a specific final state s^s^p^ can be obtained by 
computing the imaginary part of the forward-scattering amplitude using unitarity cuts 
and selecting only cuts corresponding to the desired final state |3lj . At leading order, this 
simply amounts to multiplying by the branching ratios BR(gx — > q[pi) and BR(g 2 — > ^2)- 
The matching coefficients C pro d of the production operators are computed from the on- 
shell amplitude for the processes pp' —> 5i§2 at the desired accuracy, as sketched for LO in 



2 As written, the forward-scattering amplitude is not infrared safe but requires mass factorization where 
matrix elements of the collinear fields arc identified with parton distribution functions (PDFs). 
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Figure [3j The coefficients C± p are calculated by expanding the forward-scattering amplitude 
in the hard momentum region, where all loop momenta are far off-shell]^] The leading non- 



resonant contribution to the processes (2.5) arises from the squared tree amplitudes of 



the processes pp' — > §is' 2 P2 and pp' — > s[pis 2 as sketched in Figure [4] for the example of 
squark-antisquark production and decay. The EFT calculation of the forward-scattering 
amplitude to order a s and 5 includes soft one-loop corrections, the one-loop corrections to 
Cprod as well as kinematic corrections due to subleading kinetic Lagrangian terms. 

2.3 Soft and Coulomb resummation for the resonant contribu- 
tions 
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The resummation of soft logarithms and Coulomb corrections has been derived in 
by establishing a factorization of the doubly resonant contribution to the partonic cross 
section into hard, soft and non-relativistic (potential) matrix elements, 

= 52H%(m 9 ,m 9 ,n) J du J Ra (E + if - ^)W R -(u,fx) . (2.12) 
Here we have defined the average width 

f = ^(r 5l + r s - 2 ) (2.13) 

and R a are the irreducible colour representations in the decomposition of the product of the 
colour representations of the sparticles s\ and §2- In analogy to the zero-width treatment 
in [12], as a default we use the expression 

E = fh(l- 4 ^) (2.14) 



for the energy of the squark or gluino pair which coincides with the non-relativistic expres- 
sion E = fh/3 2 for s > Am 2 . Therefore we refer to this treatment as the '/3-implementation'. 
In order to estimate ambiguities of the threshold approximation we will also use an l E- 
implementation' defined by E = y§ — 2m that agrees with ( 2.14[ ) near threshold. 



The hard functions H^t are related to the square of the matching coefficients C pro d. At 
LO they are proportional to the Born cross section at threshold]^] The soft functions W Ra 
are matrix elements of soft-gluon Wilson lines and implement the eikonal approximation. 



Precise definitions can be found in 11 . Threshold resummation of soft logarithms can be 
performed using renormalization group equations for the hard and soft function [9,54,55 
that can be solved in Mellin-moment or momentum space 56 . The relevant anomalous 



■52 



53 



3 Explicit calculations have been performed for W^-pair production 31 and top-pair production 
at electron-positron colliders. 

4 As in [12], in our numerical results we compute the hard function from the full Born cross section 
instead of its threshold limit. 
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dimensions for all squark and gluino production processes up to NNLL accuracy have been 
collected in (9 11 



The potential function J Ra is an expectation value of the non-relativistic fields ip(°\ 
that have been decoupled from the soft gluons by a field redefinition and interact only 
through the exchange of Coulomb gluons. It is related to the zero-distance Green function 
of the non-relativistic Schrodinger equation with a Coulomb potential, G^ 1 (0,0; £): 



J Ra (£) = 2imG£ QW (0,0;£) 



R a (0), 



(2.15) 



with the Coulomb Green function including all-order gluon exchange in the MS-scheme 
58l: 



57 



G* a{0) (0,0;£) 



(2m red )' 

47T 



£ 



2m red 



[-D Ra )a s 



1 



In 



m red £ \ 
V 2 J 



-^+7^ + ^(1- 



2y/-£/(2m red )J 



(2.16) 



Here je is the Euler-Mascheroni constant and m red = fn^m^/^m^ +m s - 2 ). The coefficients 
D Ra = §(Cr q — C Rl — Cr 2 ) of the Coulomb potential depend on the quadratic Casimir 
operators of the colour representations Ri of the sparticles Sj and of the irreducible repre- 
sentation R a in the decomposition of -Ri <g) R^- The trivial (a s — > 0) potential function for 
finite width obtained by neglecting the Coulomb corrections is given by 



v^^ + f 2 ) 1 



e + (e 2 + r 2 ) 1 / 2 



TT 



(E+ (E 2 + f 2 ) 1 /2) 2 + f 2 



(2.17) 



The use of a complex energy E + iT takes the finite-width effects in the doubly-resonant 
cross section into account in leading order in the non-relativistic expansion [17]. 

For energies below the production threshold E < and for F = 0, the Coulomb Green 
function develops a series of bound-state poles at energies 



E n 



2m red a 2 s D 2 Ro 
An 2 



(2.18) 



These contributions have been taken into account in the NLL predictions of |12j. For 



unstable particles, the complex energy shift in the potential function in (2.12) leads to a 



smoothing of the discrete set of bound-state poles (2.18) into a continuous non- vanishing 
contribution below the nominal production threshold. 



3 Estimate of contributions to the cross section 



In Section [2] we have discussed the treatment of finite- width effects (Section |2.2 ) and 
higher-order soft and Coulomb corrections for the case of stable particles (Section |2.3). 
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The resummation of threshold logarithms and Coulomb corrections was performed under 



the assumption (1.2) that a sizable contribution to the total hadronic cross section arises 
from the region where a s ln/3 ~ %■ ~ 1, leading to a reorganization of the perturbative 



expansion according to (1.3). For finite decay widths, the replacement E — > E + iT in the 



potential function leads to the screening of the threshold singularities, raising the question if 



the resummation based on the counting (1.2) is still justified. Furthermore, for an unstable 



particle the size of the non- resonant contribution to the cross section (2.7) compared to the 



higher-order corrections to the resonant contributions has to be addressed. In this section 
we use power-counting arguments to estimate these effects and compare these expectations 
to explicit results for squark and gluino production. 



3.1 Size of the decay width 

To discuss the size of the screening effect and the non-resonant contributions, we consider 
the decay of the sparticle s of mass M into a sparticle s' of mass m and a massless parton 
p, mediated by the strong interaction. On kinematic grounds, the decay width is of the 
form 

r(S ->■ s'p) = a s M (1 - x 2 ) f{x) (3.1) 
for some function / of the mass ratio 



x 



m 
M' 



(3.2) 



In the following, we consider the decay width (3.1) for two limiting cases: 
a) Light decay product, x — > 0: 

V/M ~ a s . 



(3.3) 



In this case, the expansion parameter S of the effective theory (2.9) is estimated as 



<5~/3 2 



a.. 



(3.4) 



b) Heavy decay product, x — > 1: 



T/M ~ a s x (1 - x 



2N7 



(3.5) 



where 7 = 1 if the function f(x) in ( |3.1 ) is of the order one, but 7 = 2 for the case of 
squark and gluino decay (2.4). In this case, another small scale k = (1 — x 2 ) 7 is present, 



and the hierarchy of scales is given by 

5 ~ (3 2 ~ a s n. 



(3.6) 



In particular, if for squarks and gluinos the numerical size of the mass ratio is such that 
K = (1 - (m/M) 2 ) 2 ~ 0.1, i.e. m/M > 0.8, we count 5 ~ a 2 . 
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For the phenomenology of unstable strongly-interacting particles, different scenarios 



have been identified (e.g. 21 59]). For very narrow particles, bound states can form that 
subsequently decay into di-photon or di-jet final states and therefore lead to a very different 
phenomenology than the missing energy signatures usually employed in squark and gluino 
searches. In the MSSM this case is only relevant for gluinonium and stoponium production 
and will not be considered further here. For a sparticle decay width that is larger than the 
bound-state decay widths into di-photons or di-jets, but much smaller than the energy E\ of 



the first bound-state (2.18), a few bound state-peaks can form, but the decay is dominated 
by the underlying sparticle decay, while for T < Ei, bound-state effects enhance the cross 
section near threshold, but the different peaks might not be separated. Since Ei ~ a 2 s , 
bound-state peaks might be visible in the invariant mass-distribution below threshold for 



narrow sparticles in the category b), as we will study in Section 4.1 For Y > E\ the 
bound-state effects are washed out, which we expect to be the case for broader sparticles 
in category a). 

We will now estimate the size of the Coulomb and soft corrections as well as the non- 
resonant contributions for these two limiting cases. 

3.2 Coulomb corrections 

In the two scenarios for the decay width, the n-th Coulomb correction is of the order 

\ n / 2\ n/2 ( rn/2 ™/ 2 \ 

a s\ I a s\ I ° as , case a) 



PJ \5J \{f) n/ \ caseb) (3 ' ?) 

Therefore by power counting, Coulomb resummation is not necessary for case a), i.e. for 
particles with decay widths of the order of Y/M ~ a s ~ 10%. This exemplifies the screening 
of the threshold corrections by the finite lifetime. However, to reach NLO accuracy in 5 ~ 
a s the second Coulomb correction has to be included. This is analogous to the case of Im- 
pair production considered in |3Tj (up to the replacement of a by a s ). In the opposite limit 



of narrow particles where numerically ~ o? s ~ 1% (case b), k ~ a s so the screening of 
the Coulomb corrections is not effective and Coulomb resummation should be performed as 
in the stable casej^] Since in these formal counting arguments the numerical prefactors such 
as the strength of the Coulomb potential are not taken into account, the actual relevance 
of the Coulomb corrections has to be studied on a process-by-process basis. In Figure [5] 
we show the contribution of the first (Ci, red dotted line) and second (C*2, blue dashed 
line) Coulomb corrections, as well as the resummed Coulomb corrections (Cf u n, solid black 
line), to the total hadronic production cross sections as a function of the decay width for 
all four squark and gluino production processes. The Coulomb corrections are normalized 
by the approximation <x^ resh that is defined as the effective theory cross section (2.12) 



without soft-gluon resummation and using only the a° term of the potential function (2.17) 



One sees that for larger widths the Coulomb corrections are increasingly saturated by the 



5 This is analogous to the well-known case of top-quark production at linear colliders, where the small 
decay width of the top is due to the electroweak nature of the decay, T t /m t ~ chew ~ 
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PP -> q q , r = 0.5 



Ao "c / o-n,,„i, 





f[ GeV] 




r[GeV] 0.00 1 



r[GeV] 



Figure 5: Size of the first Coulomb correction (red, dotted), first and second Coulomb 
corrections (blue, dashed) and resummed Coulomb corrections (black, solid) relative to 
the leading cross section in the effective theory as a function of T and for fh = 1500 GeV. 
r is defined as r = mg/rriq. For technical reasons the potential function is set to zero for 
E < AE as defined in p~6j). 



first two Coulomb corrections, while the second Coulomb correction always gives a non- 
negligible contribution, in agreement with the general discussion above. However, the total 
contribution of the Coulomb corrections can be sizable, as for gluino-pair production, or 
numerically small, as for squark-squark production. 

3.3 Soft corrections 

The parametric scaling of soft corrections in the scenarios a) and b) considered in Section 



3.1 is given by 



ln/3~ln^~ \\ lna - CaS6aJ (3.* 

I ^[lna s + In k) ~ lna s , case b) 

where in case b) we have again considered the case where numerically k ~ a s . It is thus 
not obvious at a first glance in which scenarios resummation of soft-gluon correction is 
necessary, if at all, since for a s ~ 0.1 one has a s \ lna s \ ~ 0.23. Clearly one would expect 
resummation effects to be more important for case b), where the screening effect due to 
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PP -» q q , r = 0.5 PP -» q q , r = 0.5 




Figure 6: Size of the NLL soft corrections expanded to NLO (red, dotted), up to NNLO 
(blue, dashed) and the resummed NLL corrections (black, solid) relative to the leading cross 
section in the effective theory as a function of T and for fh = 1500 GeV. The remaining 
definitions are as in Figure |5j 



a finite width is less effective. However the parametric scaling in the two scenarios is 
formally the same and equal, up to constant prefactors, to ~ lna s . As already stressed 
for the case of Coulomb corrections, the scaling argument which leads to (3.8) does not 
take into account the coefficients of the logs, which can be numerically large. This is 
particularly true for gluon-initiated squark-antisquark production and for squark-gluino 
and gluino-gluino production, where the Casimir invariants for the adjoint and higher 
colour representations appear in the prefactors. Once again, a meaningful assessment 
of the importance of resummation can only be done on a process-by-process basis by a 
numerical analysis of different fixed-order contributions. 

Figure [6] shows the effect of NLL soft resummation (Sf u u, solid black line) on the cross 
section of the four SUSY-pro duct ion processes considered in this work and the contributions 
of the 0(a s ) (Si, red dotted line) and 0(a 2 s ) (£2, blue dashed line) terms obtained from the 
expansion of the resummed cross section, as a function of the width f . As in Figure [5] the 
cross sections are normalized to the leading threshold approximation. It can be seen that 
for small widths of order T/m ~ a 2 s resummation is numerically important for all processes 
(and more so for those where large colour charges are involved), with NNLO corrections of 
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the order of 10 — 25% of the LO cross section and terms beyond NNLO as large as ~ 10% 
for squark-gluino and gluino-gluino production. One can therefore argue that in scenario 
b) it makes sense to keep the zero-width scaling a s In (3 ~ 1 for soft logarithms, which leads 



to the representation of the cross section given in (1.3). 

At larger width (T/rh > a s ) the bulk of NLL corrections is accounted for by the 0(a s ) 
terms, which correspond to a 30 — 50% contribution to the Born result. However NNLO 
terms still represent a correction of order 5 — 15% depending on the process considered. 
Consistent with the stronger screening effect expected for larger widths, higher-order soft 
corrections beyond NNLO are small, < 5% for gluino-gluino and squark-gluino production 
and below 1% for squark-antisquark and squark-squark production. We thus conclude that 
in case a) the relevance of soft resummation is not immediately clear, though the inclusion 
of higher-order corrections at NNLO might still be necessary to achieve an accuracy of 
order 5 ~ a s . It is interesting to note that for large widths the fixed-order one-loop soft 
correction is numerically about ^Ja~ s ~ 30% of the tree- level threshold result, and the 
two-loop soft terms of order {y/a~ s ) 2 ~ 10%. 

To summarize, in the two limiting scenarios one can give the following parametric 
representation of the cross section: 

a) Adopting for the soft logarithms the (approximate) scaling a s In 2 j3 ~ y/5 ~ y/oT s , which 



is motivated by the numerical results in Figure |6j and using (3.7), the cross section can 
be represented as 



n 



<5pp/ oc x < 



a. 



a s In /3 

2 

J 



LO 

N^LO 

{a s In 2 /3) 2 , (a s In (3) 2 /(3, a s ln/3, /3 2 NLO 



(3.9) 



where the expansion in half-integer powers of 5 is similar to the case of W-pair produc- 
tion at a linear collider 31 . There are no terms linear in (3 because they are known to 



average out to zero for the total cross section. 



b) For the case k 



the counting is identical to the stable case. 



While according to (3.9) resummation is parametrically not necessary for T/rh > a s , 
in practice it is not always clear for which numerical value of T one can switch from 
the prescription (1.3) to (3.9). This is particularly true for the transition region where 
a 2 < T/rh < a s . For this reason in Section we will use the NLL implementation of (1.3) 



for arbitrary values of the width. This choice has the advantage of including, in both limits 
of a small and large width, all the terms relevant to achieving an accuracy of ~ 5, with 
the exception of the /3 2 contribution appearing at NLO in the counting (3.9). These are 



however correctly taken into account by matching the effective-theory result to the exact 
Born result including the leading finite-width and power suppressed effects, as explained 
in Section |3]4l 
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3.4 Non-resonant corrections 



We now turn to an estimate of the size of non-resonant corrections neglected in present 
higher-order calculations that treat squark and gluinos as stable. The leading non-resonant 
contributions as well as subleading kinetic corrections to the double-resonant cross section 



are both included in a calculation of the full Born cross section of a process of the form ( 2.5 ) 



using multi-leg Monte Carlo programs [24] so the strict EFT treatment is not necessary. 
Also the non-resonant diagrams are in general well-defined only with kinematic cuts, as 
mentioned above. A study of the full non-resonant effects for all squark and gluino pro- 
duction processes would therefore require a dedicated study taking the realistic selection 
cuts used by the LHC experiments into account, which is beyond the scope of this work. 
We will therefore limit ourselves to an estimate of the unknown radiative corrections to 
the non-resonant contributions that cannot easily be computed in a Monte Carlo program. 
These are given by higher-order corrections to the diagrams shown in Figure |4j The com- 
putation of these corrections is currently an ongoing effort for top-pair production at a 



linear collider 52 53 60 61 and the computation for the case of squarks and gluinos is 
far beyond the scope of the present work. These corrections would also be included in a 



complete fixed-order NLO calculation of the process (2.5) 



We consider again the two cases introduced in Section 3.1 



a) Since in the calculation of the matching coefficients of the four-parton operator all the 
momenta are large compared to the scale S, there is no enhancement of these corrections 



by resonant propagators and one expects them to be of the order 31 



^non-res _ 3 i.res - , 
pp s w pp' 



1 



(3.10) 



For the scaling (3.4) appropriate for T/M ~ 10%, the leading non-resonant corrections 
are therefore of the order ~ yjY~/M, as the first Coulomb correction. The unknown 



radiative corrections to the non-resonant contributions (3.10) therefore are of the order 



a 



~ a 



3/2 



relative to the leading resonant cross section, and therefore beyond 



NLO accuracy. 



b) For a small mass hierarchy leading to the scaling (3.6) the non-resonant contributions 
can be further expanded according to /3 ^J~p <C 1, with p = 1 — m/M. For the case 
of top-quark production the leading non-resonant corrections were found to be of the 
order 



53 



Using (3.5) we have 



r 



"iion-rcs "res . . 
°pp' ~ a pp' X 



a s (l-x 2 ) 7 



Mp 



(3.11) 



(3.12) 



Mp (1-x) 

For squark and gluino decay 7 = 2, so counting n ~ a s we find that the suppression of 
the non-resonant corrections is of the order of (tj— rv>) 1 ^ 2 ~ (T/M) 3 / 8 compared to the 



• Mi/2, 



14 




Figure 7: The partonic cross section for squark-gluino production and squark decay as a 
function of the energy E = m(3 2 at a fixed average mass fh = \{rriq + nig) = 1500 GeV. 
Left: r = ^ = 0.95, f = 0.4 GeV. Right: r = 0.5, f = 32 GeV. 



leading resonant term. The unknown radiative non-resonant corrections are then of the 
order (^a^ 2 ) 1 / 2 ~ (r/M) 7//s , i.e. of a similar magnitude as for the case x — > 0. These 
corrections are beyond NNLL accuracy. 

To study the numerical impact of the non-resonant corrections, we consider the following 
approximations: 

• ^nwa : The LO cross section calculated using the narrow-width approximation ev- 
erywhere, i.e. the on-shell production cross section in the r — > limit is multiplied 
by the decay branching ratios. This approximation is valid for M ^> T and far above 



threshold [25]. In this work only the decay processes (2.3) are considered and thus 
4wa = <%S r = 0)- 

^Thresh- The threshold approximated LO cross section in the EFT, using only the a° 



term of the potential function (2.17). The unstable-particle momentum in the centre 



of mass frame is in the potential region (slightly off-shell). 

(Tq^: The LO cross section obtained using the off-shell doubly-resonant diagrams and 
fixed-width Breit-Wigner propagators computed with WHIZARD [62] and validated 
with SHERPA |63| when possible. This implements tree- level off-shell effects and 
allows for unrestricted momenta of the unstable particles^ 



In Figure [7] the three different approximations of the partonic cross sections are shown for 
the example of squark-gluino production and subsequent squark decay as a function of the 
energy E = faff 1 for two examples of the decay width. The plots for the other processes 



6 Thc selection of doubly-resonant diagrams violates gauge invariance. However, we have verified for 
a selection of different covariant and axial gauges that gauge violation is numerically below 1% for the 
processes of interest in this work. 
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have the same behaviour as above and are not shown here. It can be seen that the EFT 
approximation ^Thresh a g rees with the off-shell results for small widths, but is systematically 
shifted for larger decay width. This is the expected effect of the non-resonant contribution 
to the Born cross section, that is of the form <r°°, n " res ~ 



m 



31,52,60 



pp 



ImC 4p , and has also been observed 



Note that while we do not include truly singly or non-resonant diagrams 
in the off-shell approximation, the off-shell momentum configurations of doubly-resonant 
diagrams contribute to the non-resonant part of the cross section and, in fact, provide the 
dominant numerical effect in the case of W or top production 31 52 53 . Therefore the 



shift observed in Figure |7| provides an estimate of the non-resonant corrections. For T/rh = 
32/1500 = 0.02 the nonresonant correction is of the order of 40%, which is of the same 
order of magnitude, but somewhat larger, than the estimate (r/M) 1 / 2 -(r/M) 3 / 8 = 14%- 
24%. We then estimate that the uncalculated corrections due to higher-order non- resonant 
effects are of the order of < 5%. 



4 Application to squark and gluino production at LHC 

In this section we present results for finite-width effects on the soft and Coulomb NLL 
corrections to squark and gluino production at the LHC in the setup discussed in Section [2] 



In Section 4.1 we give results for the invariant mass distributions, while in Section |4~2| we 



discuss our implementation and present numerical results for the total cross sections. 
4.1 Invariant-mass distributions 

In order to discuss the finite- width effects on the production cross sections, we first consider 



the invariant-mass distribution. Similarly to the total cross section (2.12), it satisfies a 



factorization formula near the production threshold [IT] (see also e.g. |21| " 

=^^(M glS - 2 -2m+.r) 5>i^(/i) [ 1 dTL ppl (r^)W R -(2(V~s-M SlS2 )^), 

SlS2 Rc P,p'=q,q,9 Jto 

(4.1) 

where r = Mj ' s Is. We have defined the parton luminosity in terms of the PDFs f p /N for 
a parton p in a nucleon N: 

Lpp>(T,n)= dx x dx- i S(x 1 X2-T)f p / Nl (x 1 ,ii)f p/ / N2 (x2,^). (4.2) 



In our numerical results we use the NLL-resummed soft function that can be found in 12 
Since we are interested in the invariant-mass spectrum near threshold we deviate somewhat 
from the default treatment of the total cross section and always express the hard function in 
terms of the threshold-approximation of the Born cross section and use a constant Coulomb 
scale, defined as solution to the equation fi c = 2a s (fic)m r \D Ra |, and constant soft scales 
Us that are process-specific and determined in (12). 
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Consistent with the treatment of the total cross section discussed in Section 12.31 as a 



default we use the expression (2.14) for the energy of the squark or gluino pair, i.e. we 



write the argument of the soft function W Ra in (4.1) as 



(4m 2 \ 
1 - — J + 2m-M Sl ~ S2 . (4.3) 

As for the total cross section, this will be referred to as the '/?- implementation'. Note that 
the soft function vanishes for negative argument which implies that the lower boundary of 
the r integral is given by 

4m 2 m 
s dm — M Sl g 2 

In order to estimate ambiguities of our treatment, we also introduce an '^-implementation' 
keeping the unexpanded expression Vs — Mg 1 s 2 in the argument of the soft function. In 
order not to introduce artificial differences in the two implementations due to a different 
boundary of the r integral, we perform the replacement Mg 1 s 2 — > M~ s = m(3— 4m 2 /M| i ; 
in the ^-implementation and consider the differential cross section with respect to M' Sii 

Based on experience with VT-pair production [31] we do not expect the leading resonant 
approximation to be valid significantly below the nominal threshold, where higher-order 
terms in the non-relativistic expansion and non- resonant contributions become important. 
The differences between the default (5 and the E approximation provide an estimate for 
the size of the higher-order relativistic effects. Note that because of soft radiation (as 
taken into account by the convolution of the soft function with the parton luminosity) 
also the invariant mass distribution directly at threshold is sensitive to larger partonic 
centre-of-mass energies and the associated ambiguities in the threshold expansion. 

In Figure [8] the invariant-mass distributions for the four relevant squark and gluino 



1S2" 



production processes (2.6) are shown for m = 1500 GeV at a centre-of-mass energy y/s 



TeV for two different mass ratios corresponding to the ) and b) discussed in Section 



3.1 . For the case m- g > rriq we consider r = mg/rriq = 2 and r = 1.05, corresponding 
to Yg/rrig = 12% and Yg/rrig = 0.2%, respectively. For the case rrig < rriq we consider 
r = 0.5 and r = 0.95, corresponding to Yq/rriq = 3% and Yq/rriq = 0.05%, respectively. As 
explained in Section |3j we will perform treshold resummation also for the case of larger 
decay widths, although it might not be strictly necessary in this case. 

In the small-width case, i.e. the left-hand plots in Figure [8] where r = 0.95 or r = 
1.05, the ambiguities in our prediction as estimated by the difference of the E- and f3- 
implementations are moderate. Depending on the process, the different regimes mentioned 



in Section 3A_ are observed: for the given choice of parameters, a single would-be bound- 
state peak appears for squark-squark production, where f = 0.9 GeV < E± = 1.5 GeV, 
while two peaks are visible in gluino-pair production where T = 3 GeV <C Ei = 30 GeV. 
In the remaining two processes the second bound-state peak is less pronounced, but still 
visible. 

In the case of larger widths, i.e. the right-hand plots in Figure [8] with r = 0.5 or r = 2, 
the decay widths are always much larger than the bound-state energy, so no peaks below 
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Figure 8: Invariant- mass distributions for an average mass fh = 1500 GeV of the produced 
sparticle pair and a centre-of-mass energy y/s = 8 TeV for the (3 (blue, solid) and E (red, 
dashed) implementation as defined below ( 4.3[ ). Left-hand plots are examples of case b) 
{Y /fh ~ right-hand plots case a) {Y /fh ~ a s ). 
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threshold are observed. It can be seen that the invariant-mass distributions do not vanish 
for M Sl s 2 far below the nominal threshold, M SlS2 = 3000 GeV. This behaviour arises as a 
combination of the rise of the PDFs for small x and the behaviour of the potential function 



m 



4.1 



that scales as J Dr oc T/(2m - M Sl S2 ) 1/2 ->■ r/(2m) 1 / 2 for M hh -> 0, which follows 
from the LO (trivial) potential function (2.17). While for the default /3-implementation 
the distributions approach a flat plateau far below threshold, in the ^-implementation one 
observes a (clearly unphysical) rise for larger widths for processes involving gluinosQ In 
fact, as already pointed out, far below threshold the LO potential function is not expected 
to be a good approximation, and the inclusion of non-resonant contributions and relativistic 
corrections are needed for a realistic description of the invariant mass distributions in this 
region, similar to what was seen for the partonic cross sections in Figure [7} 

The difference between the E- and /^-implementations is amplified due to the soft 
function in (4.1), W Ra (u) ~ u~ 1+2v with rj < 0. In the ^-implementation the argument of 
the soft function (left-hand side of Eq. (4.3 )) approaches zero faster than the corresponding 
quantity in the /^-implementation (right-hand side of Eq. (4.3)) as the lower integration 
boundary ro is approached, and the soft function W diverges faster]^] If soft corrections 
are not taken into account, the soft function reduces to a delta function, W Ra {oj) oc 5{oj). 
In this case we have observed a substantial reduction in the difference between the E- and 
/3-implementation also for large widths. 

From the results presented in this section we conclude that our approximation is rea- 
sonably accurate for processes with decay widths of the order of T/rh < 5%. As we will 
argue in Section 4^2, this covers the phenomenologically-relevant MSSM parameter-space 
regions. Our effective-theory treatment becomes inadequate for larger decay widths, so 
our results in this regime should be considered as indicative only. 



4.2 Total cross sections 

In order to compute the total cross section for a finite decay width, in principle the u- 
integral in (2.12) has to be performed up to the maximal energy oj = 2\/§ since the 



potential function J^ a is defined for arbitrary negative values of the partonic energy E. 
However, as mentioned above, far below threshold (i.e. u ^> 2E) the leading threshold 
approximation becomes insufficient and higher-order and non-resonant contributions be- 
come relevant. This was shown by the difference of the o^hLsh and ctq^ curves for large 
widths in Figure [7] and by the large difference of the /3- and ^-implementations for the 
invariant-mass spectrum and the unphysical behaviour of the latter below threshold, as 
seen in Figure |8j Since the non-resonant and higher-order non-relativistic corrections for 
squark and gluino production are not available beyond LO, we introduce a cut-off in the 



7 Note that this is not an artifact of our use of M' s s for the ^-implementation. In fact, for the original 



formula (4.1) an even stronger rise is observed, since smaller x values are probed in the parton luminosity 
due to smaller values of To- 

8 This is further enhanced by the rise of the PDFs at the lower integration boundary. 
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potential function: 



J R {E + iT) = 2lm[G R (E + tT)]6{E + AE) , AE > 0, (4.5) 

with E equal to the expression ( |2.14[ ) as our default implementation. The parameter AE 
cuts off the convolution of the soft and potential function at u — 2(E + AE). In order to 
capture the bound-state peaks, this cut-off is chosen as 



AE = l.l|£a| + AE, (4.6) 



where E\ is the energy of the lowest-lying bound-state peak, see Eq. (2.18). The prefactor 



1.1 multiplying E\ insures that the bound-state contribution is always fully included. As 
a default for AE we take 



AE := MinflOT, 0.25m]. (4.7) 



A related cutoff AE = \Ei\ +3T was used for gluino-pair production in 21 where the focus 
was on narrow gluinos with decay widths at most of the order of Tg ~ \E\\ ~ m g a 2 s . In 
the case of larger decay widths, the value of AE = 10f becomes of the order of m§, where 
we expect the threshold approximation to become invalid. Therefore we have chosen the 
cutoff AE ~ 0.25m in the case of larger decay widths. The uncertainty in choosing AE is 
estimated by varying it around its default value: 



AE /2 < AE < 2AE . (4.8) 
The hadronic cross sections are obtained by convoluting the partonic cross sections 



with the parton luminosity (4.2) 



VniN 2 ^ss'x{s) = j dr ^ L pp'( t i f J 'f)vp P '( TS i Vf) > ( 4 - 9 ) 
T ° p,p'=q,q,9 

where the lower integration boundary depends on the cutoff, 

/ » \ 4m 2 1 

r AE = rrr^- 4.10 

uv ; s 1 + AE/m v ; 

In the ^-implementation, the lower integration boundary is instead given by t (AE) = 
(2m — AE) 2 /s. In order not to introduce numerical differences due to different integration 
boundaries, in the ^-implementation we use the same numerical value for To as for the 
/3-implementation, and adjust the value of AE in the potential function accordingly^] 

While the Born amplitudes for the production and decay processes of squarks and 
gluinos can be treated exactly with Monte Carlo programs [24], the NLO corrections for the 
production processes are known only in the narrow-width approximation, as implemented 
in Prospino |64|. In order to take these exactly-known contributions into account, we 



9 This is related to the replacement of the invariant mass Mg 1 s 2 —> M'~ - for the ^-implementation in 



Section 4.1 
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match our finite-width NLL results to the sum of the Born cross section including finite- 
width effects and the NLO corrections in the narrow-width appr oxim ation. In practice, 



we approximate the Born cross section by CTq^ defined in Section 3.4 and computed with 
Whizard. The total hadronic cross sections at NLO for finite widths are thus approximated 
as: 

0NLo(r) = CTofftT ) + (ct NLO ,NWA - 4wa)' I 4 - 11 ) 

Here ctnlo.nwa is the NLO result obtained with Prospino and ct^wa i s the leading-order 
production cross section for stable squarks or gluinos (see Section 3.4). Analogously to the 
narrow-width case [l2], the NLL resummed cross section is then matched to CTNLo(r): 

^NLL hcd (f) = Aa NLL (r) + a NLO (r), 

Acr NLL (r) = a NLL (r) - 0- NLLj LO (r) - ct NLLi nlo (r = o). (4.12) 

Here ctnll, (n)lo is the expansion of the NLL resummed cross section to (N)LO accuracy. 
For comparison we will also consider results where the decay width is set to zero in the 
higher-order corrections, i.e. 

a matched (0) = A^Q) + ^(f ) . (4.13) 



The result (4.12 ) is our best current prediction. It could be made more realistic by replacing 
the off-shell approximation for the tree cross section by a full LO simulation including 
acceptance cuts, which we leave for further investigation. Contrary to the NLO cross 
section o"NLo(r), the matched NLL cross section ctnll^ ) does contain finite-width effects 
at NLO from the resummed Coulomb and soft corrections. It is also interesting to consider 
only the finite-width corrections beyond NLO. This shows the size of the effects that are 



beyond a possible future exact NLO calculation of the processes (2.6). For this purpose 
we consider the alternative matching 

^NLL ChCd ' (2) (f) = AaS L (f) + a NLO (r), 

Acr^£ L (r) = ctnll (r) - ctnll, lo(F) - o"nll, nlo (T) . (4.14) 

For the results shown in this section we take as a representative average mass m = 
1500 GeV and the setting is the LHC at a centre-of-mass energy of y/s = 8 TeV. We 
use the /3-implementation as default and the scales are taken to be the same as in the 



zero- width case 12 : fif = m = 1500 GeV, fih = 2m = 3000 GeV and a running soft and 



Coulomb scale. We note that the soft and Coulomb scales are always fixed below threshold. 



The LO widths in Eq. (2.4) for the squarks and gluinos are used. 



Figure M shows the ratio of the matched NLL cross section for finite widths (4.12) 



to that for the zero- width case (4.13) (red, dot-dashed) as a function of the gluino-to- 
squark mass ratio. The green band gives the total uncertainty for the zero-width NLL 
result, including ambiguities in soft, hard, factorization and Coulomb-scale choices and the 
difference of the E- and /3-implementation, added up in quadrature as discussed in detail 
in 12 . The difference of the E- and /3-implementation alone is given by the blue band. 
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Figure 9: The plots show the ratio cr^ ched (r)/cr^ ched (0) of the matched finite-width 



NLL cross section (4.12) to the zero-width approximation (4.13) (red dot-dashed, left- 
hand vertical axis) for m = 1500 GeV and a/s = 8 TeV and the ratio T/m in per-cent 
(purple dotted, right-hand vertical axis). Both quantities are plotted as a function of the 
gluino-to-squark mass ratio. The NLL result is for default AE = AE . The red band 
represents the error from AE variation and E vs (3 uncertainty added in quadrature, the 
green band the total error for zero width, the blue band the E vs uncertainty for the 
zero-width case. 



The red band indicates the ambiguities related to our impleme ntat ion of the finite-width 
effects, as estimated by the AE variation as described in Eq. (4.8) and the difference of 
the E- and /^-implementations, summed in quadrature. The plots also show the ratio of 
the average decay width to the average mass of the produced particles. 

It can be seen that for all processes the finite-width result including uncertainties lies 
within the uncertainty estimate of the NLL calculation for stable squarks and gluinos if 
f/m < 5%. For larger widths, large corrections are observed for the gluino-gluino and 
squark-gluino processes, where the ratio T/m grows up to 12%, and the finite-width cor- 
rections become of the order of 20%. In this case the uncertainties due to the treatment of 
finite-width effects are of a similar magnitude as those due to scale and resummation ambi- 
guities. This indicates that corrections beyond the leading E — y E + iT replacement, such 
as higher-order non-relativistic corrections and non-resonant corrections become relevant. 
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Figure 10: Same as in Figure 9j but for the ratio cr^^ cd) ^ (V) / a^^ hed {Q) of the matched 



cross section (4.14) including only finite-width effects beyond NLO to the zero-width 



approximation (4.13). See Figure 9 and the text for explanation. 



A more accurate treatment would require the matching of the NLL results to the exact 
NLO QCD cross sections for the full qqq or gggg-processes and the inclusion of the N 3 / 2 LO 
corrections in the counting (3.9) that are expected to be of the order <5 3//2 ~ (f /m) 3 / 2 ~ 5%, 
both of which are not available at the moment. Nevertheless it is interesting to anticipate 
such a future matching to a full NLO calculation by using the matching prescription (4.14) 
where the finite- width effects are included purely beyond NLO. The results are shown in 
Figure 10, where the colour coding used for the lines is the same as in Figure [9j It is 
seen that the finite-width effects on the central values, as well as ambiguities due to the 
AE variation and the E- and /3-implementation difference, are much smaller than for the 
matching (4.12). In particular, the total uncertainty bands related to the finite-width 
implementation (red) lie now fully within the total zero- width NLL error bands (green). 
The difference of the plots in Figure [9] and Figure 10 indicates the potential impact of a 
full NLO SQCD calculation for the processes (2.6) and shows that the finite-width effects 
beyond NLO are expected to be small. 

To assess the size of the finite width corrections relative to the NLL soft and Coulomb 



corrections, in Figure 11 we compare the K-factors for our NLL result including finite de- 
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Figure 11: The i^-factors of the NLL cross sections including (red dot-dashed) and ex- 
cluding (black) the widths for default AE = AE and rh = 1500 GeV at centre of mass 
energy y/s = 8 TeV. Dashed black: NLL resummed cross sections excluding bound-state 
contributions. Green band: Total error for zero width. 



cay widths to the result for stable squarks and gluinos (4.13). The i^-factors relative to 
the NLO prediction are defined as -?Tnll = onll/ctnlo, where the finite- width approxima- 
tion (4.11 ) is always used for the denominator. We also show the result NLL no Bs defined as 
the cross section in the stable case, excluding the bound-state contribution below thresh- 
old. In general, for relatively narrow particles with T/m < 5% the K-factors including 
finite-width effects are well approximated by those for the zero-width case including the 
bound-state contribution, with the exception of squark-gluino production and gluino-gluino 
production for the case rrig > rriq, where larger finite-width effects are observed. 

The relatively large finite-width corrections for gluino-production processes for rrig > rriq 
have to be compared to the contributions of these processes to the total SUSY production 
rate. This can be seen in Figure 12, where in the left plot the i^NLL-f actor is shown for 
the finite-width and zero- width case. The effect of the QCD widths is seen to be negligible 
and the finite width can be safely set to zero when it comes to the total SUSY cross 
section. This is also expected since the numerically dominant processes are almost always 
those whose final-state particles are stable against SQCD decay, as can be seen from the 
right plot of Figure 12 The inclusion of electroweak widths is not expected to change our 
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Figure 12: Left plot: The NLL i^-factor including (dot-dashed red) and excluding (black) 
the widths for default AE = AEq and M$ = (rriq + m§)/2 = 1500 GeV at centre of mass 
energy y/s = 8 TeV. Dashed black: NLL i^-factor excluding bound-state contributions. 
Green band: Total error for zero width. The plot also shows the squark (thick orange) and 
gluino (thick green) widths (right vertical axis). Right plot: The ratio of the four relevant 
NLL cross sections to the total SUSY cross section. 
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tion of mi/ 'trig for a fixed gluino mass m g = 1000 GeV, a common light-flavour squark mass 
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Figures^] and 10 



2000 GeV. The error bands for finite- and zero-width curves are defined as in 



qualitative results. These can be dominating when r ~ 1, but are only about the order of 
T/m ~ 1% and thus the effects are expected to be negligible. 

As mentioned in Section |2j in the study of finite- width effects considered in this work 
we have neglected the decay of gluinos to stops. We do not expect the inclusion of this 
additional decay channel to significantly change the results presented in this section, except 
in the case in which gluino production is the dominant SUSY production channel and the 
gluinos only decay to stops. This is for example the case of scenarios with the mass 
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hierarchy mi < m g <C rriq, with all squarks heavy except for stops. To have an estimate 
of the effects of a finite gluino width in this particular situation in Figure 13 we show the 
NLL cross section for gluino-gluino production as a function of the ratio m^/nig for a fixed 
gluino mass m g = 1000 GeV, a common light-flavour squark mass of nig = 2000 GeV and 
top mass m t = 174 GeV. As in Figure [9] the cross section is normalized to the gluino-pair 
production cross section at zero width, and the green band represents the total uncertainty 
of the zero-width result, while the purple dashed curve gives the gluino width as a function 
of mi/ nig. As can be seen from the plot the finite- width effects amount to a maximum of 
~ 1 — 2% for small stop masses and the uncertainty associated with the finite-width result 
is contained in the zero-width case error band. We thus expect corrections to the plots in 
Figures 11, 12 to amount to at most few percents. 



5 Conclusions 

In this work we have considered corrections related to finite widths of squarks and gluinos 
on the NLL resummed results presented in |12|. As anticipated, a finite decay width leads 
to a screening of higher-order soft and Coulomb corrections. However we have found that 
for decay widths of the order of T/fh < 5%, corresponding to gluino to squark mass ratios 
of 0.5 < nig/nig < 1.4, finite-width effects on the NLL corrections are small and within the 
uncertainty of the soft and Coulomb resummation including bound-state effects performed 
in [12]. 

For larger decay widths of the order of f/m ~ a s non- resonant contributions and higher 
order non-relativistic corrections become important. In this case it becomes necessary 
to perform a calculation of order <5 3//2 ~ a s \/T /fn in unstable particle-effective theory 
matched to an exact NLO calculation for three- or four-particle final states including finite- 
width effects, which should be feasible but challenging with current NLO technology. The 
accuracy of such a calculation is expected to be better than 5%. For the doubly-resonant 
contributions to the cross section, we have observed that finite-width corrections beyond 
NLO are small. 

Because the pair production processes of the lighter coloured sparticles, which are stable 
with respect to SQCD, dominate the total SUSY cross section, finite-width effects on the 
total SUSY production rate are negligible on the whole range of squark and gluino widths 
considered. 
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